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Abstract 

We use quantum electrodynamics and the Bethe-Salpeter equation to cal¬ 
culate the bound state energies for a two-particle system comprised of a spin-0 
and spin-1/2 particle. We generalize our treatment to include the finite size 
of the nucleus and of the meson and our results can be applied to kaonic and 
pionic atoms. Our treatment includes quantum electrodynamics and recoils 
corrections to all orders. 


1 Introduction 

Exotic atoms can be of great value in the investigation of the threshold of the onset 
of strong interactions. For this it is useful to use the difference between the measured 
X-ray energy and the value of the corresponding electromagnetic energy [1]. Thus to 
gain insight into the strong interaction effects, we must remove the electromagnetic 
contribution so we need to know how to calculate this precisely. In 1990 Owen found 
a new representation of the Klein-Gordon propagator which allows a convenient for¬ 
mulation of the bound state problem using the Bethe-Salpeter equation, or which 
one of the particles is spinless 0 0 a- By separating the positive and negative en¬ 
ergy parts of the propagator, a two-component formalism was obtained in which the 
inverse of the Klein-Gordon propagator is linear in the momentum po- This makes 
it possible to combine it with the Dirac propagator, and generalize this formalism 
to include the hnite size (i.e., form factors) of both the particles. Thus our aim 
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is to calculate the energy levels of the meson-nucleus system with a purely electro¬ 
magnetic interaction which includes all recoil corrections. This necessitates the use 
of the Bethe-Salpeter equation for systems in which recoil cannot be neglected. In 
the Bethe-Salpeter equation the entire recoil contribution to the energy is included. 
For S-states the effect of the strong interaction is dominant and a comparisons with 
experiment cannot be made until this is taken into account. We leave these con¬ 
siderations for a future work since our concern here is the purely electromagnetic 
interaction. 

There has been interest in this problem due to new theoretical treatment using 
the Breit equation [5] and recent experimental results. We must point out that there 
is a fundamental difference between the Breit equation approach, which is based 
on an approximate formulation of the two body problem, and the Bethe-Salpeter 
equation which is a method for finding the location of the poles of the two-particle, 
relativistic quantum Green’s function. As it relates to quantum electrodynamics, 
this approach can give the location of the energy levels as precisely as one desires. 

2 Review of the Theory 

2.1 Lagrangian 

We take for our free Lagrangian i 0 


i - ■ -' M d^' 

£o = - rriaW + (1) 


k=l 


dx^ 2mb dx^ 


while our interaction Lagrangian is obtained from Eq. 1 by minimal coupling and M 
and jSg are given by the matrices, respectively: 


M = 


1 1 
-1 -1 


/ 3 .= 


1 0 
0 -1 


( 2 ) 


while nia and rrih are the respective masses of the particles. The total Lagrangian 
for our nucleon-meson system can be written 

£o-i/2 = Co- (3) 

where the second term is the Lagrangian of the free electro-magnetic held and third 
term represents the electromagnetic interaction between the nucleon and meson. The 
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hadronic content will be expressed by the charge distribution of these particles. That 
is, we can write [6] 


n 


0 - 1/21 


x) = -/I 0 - 1/2 = iebA'(x',t)p^(x' - x) 


2mb 2mi) 


M 


+efeA'^(x',t)p^(x' - x)<h'(a;)-—<h'(a;) - e„'0'+(x)7 ■ A'(x', t)pAr(x' -x)'0'(x) +Hc{x) (4) 

2mb 

where and are the Fourier transforms of the pion and nucleon and where the 
contribution to the Coulomb Hamiltonian is given by 





47r 


,3 ,3 t)pjv(x - xOp^(y - y') 

^ y I / /1 

|x' - y'l 


(5) 

where the integrations over x' and y' are understood and H^s describes the Coulomb 
self-interaction Hamiltonian which only plays apart when renormalization is consid¬ 
ered. In principle, we should have added an additional term to Eq. |4] to describe 
the contribution of anomalous magnetic moment of the proton to the interaction 
Hamiltonian. However, in the bound state energy levels that we consider, this term 
will not give any additional contribution so for simplicity, we have neglected it. 


2.2 Bethe-Salpeter Equation 

The Bethe-Salpeter equation [7] can be written as 

^(1, 2; 3,4) = Go(l, 2; 3,4) + Go(l, 2; 5, 6)/(5.6; 7, 8)^(7,8; 3,4) (6) 

where 

G'(l,2;3,4) = - < O|r['0'(xi)<I)'(a;2)i/''(a;3)i>'(a;4)]|O >= - < 0|T[/A(a;i)<h(a;2)'*^(x3)|)(a;4)5']|0 > 

(7) 

where 

^ dynT{nuyi) • • • 'H^nt{yn)} (s) 

The primes indicate that the helds are the fully interacting helds in the Heisenberg 
representation. Gq is the Green’s function for non-interacting helds (i.e., 4>’s and /7’s 
are not interacting with each other but the electromagnetic held is not turned oh) 
and I is the irreducible interacting kernel which is dehned by Eq. [6l Using Eq. [5] in 
the Bethe-Salpeter equation, Eq. [6] we hnd for the Coulomb kernel 
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= -i^^-io6{x^-x^)5{xQ-X8)5{xl-xl) I (fx'd^y 

47r J 1 ^ - 3 / I 

(9) 

and similarly this can be done for the other terms appearing in EqJH Except for the 
appearance of the charge densities, and Pt, these are the same as appear in the 
paper by Owen [^. So we shall not give them here. The Bethe-Salpeter eqnation 
appearing in Eq. [^can be written in C. M. coordinates 0. and for which the CM four 
vector K = {K, 0) 


,3 ,,Pn{x7 - x')p^{xs - y') 


(I)k{,x) = Gk{x, x')Ik{x', x'')(j)K{x'') where (j)K =< 0\T {i)'(pbx)^'{-pax)} \K > 

( 10 ) 

with X = Xi — X 2 and the only restriction on the p’s is that + Pb = 1- Although 
Gk{x, x') represents the free propagation of a nucleon and scalar meson in which each 
of the particles is fully ’dressed’, the order which is the accuracy of our calculation, 
we can treat each of these as free, non-interacting particles. The contributions of 
the self-interactions does not contribute to order o.^ and at most, can contribute to 
order o?. Hence we can write Gk{x,x') as0 


Gk{x,x') = 


where 


d^p 

Wv 

d^p 

w~ 


exp [-ip ■ {x - x')]SFa{VaK + p)A{pbK - p) 

^ exp [—ip ■ {x — x')] - 


2 Eb{p) {paK po)7o - 7 ■ p - rua 


■ds 


{pbK - po) - - PsTUb 


( 11 ) 




/5.= 


2 ’Tt6 




Eh{p) 2mi,Eb{p^ 


rrib+T 


ml 


\ 2mbEb{p) Eb{p) ) 


( 12 ) 


^for further details see reference [5] 

^We have modified the definition ofA(p) from what appears in [4] by a factor ds so it now reads 
A(p) = [2E{p){po — ^ — ds'm)]~^$s- This insures the correct statistics without the worrisome 
overall factor of (2po)~^ which was used in [3] 
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By expanding this in partial fractions, we can write 

Gk{x,x') = [ {x - x')]— --—Ax(x, x')^.7o (13) 

J yZn J i\ iia til) 

with 

^ n_fd^P f 1 1 1 

7(271)4 2 E{p) yr]aK+ pQ-a-p-I3ma^ r]i,K- pq-M^-( 3smb] 

(14) 

Using Eq. [TT] in the Bethe-Salpeter eqnation in CM coordinates, Eq. [101 we can 
write 

2 

[K - a--p- niaP - - mbl3s]4>K{.x) = [5(x - x') + VLk{x, x']Xloi^^d{x' - x'')6{x'q) 

2 ui}) 47r 

J Ix + u-v]^^^ ^ 

+lK{x,x')(j)K{x') (15) 

where we have written Ak{x,x') appearing in Eq{T3] 

Ak{x,x') = 6 {x — x') + ^Ik{x,x') (16) 


This will enable us to obtain a Coulomb equation that is tractable and we have 
defined 

I'j^ix, x') = Ik{x, x') - Ic(x, x') (17) 

In this work, we shall only be concerned with contributions of order and since 
it has been shown that the term VLk{x,x') contributes are at most, order so we 
can ignore this term as well. Thus, if we can treat /)^(x,x') as a perturbation to 
the zero-order equation obtained by truncating 4 from the kernel. Using our above 
treatment and the Foldy-Wouthuysen i a our zero-order equation is 


r v 2 


1 2/i j 

|x-hu-h v|_ 




(18) 


where E = K^ — rria — rrib and p is the reduced mass. 


The corrections to E we designate by AE where AE = K — Kc- To order (i.e., 

including all recoil contributions) we need to calculate both the relativistic kinetic 
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corrections (i.e., E{p) — ^ and those arising from the one photon exchange between 
the proton and the bound meson. From and from the treatment above, these are 
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AE= < 






+ 


a 


8ml Sml iml 


(TUiTVp AT (u) (v) 


(Ta • (r + u - v) X p CKTr 

+ + u 


a 


2 mamb 

1 


r + u — V 


d\d'^vpN{u)p^{v) , 
ra„ ■ (r + u - v) X p I > 


2 r + u 
r + u 


r + u — 

V 


2 ml 

-((r + u-v) ■p)p 


The first two terms in Eq. [11] are the relativistic kinetic terms, the next group of 
terms in the curly brackets are from the spin-orbit interaction and the terms in the 
last bracket are from the one-photon exchange. 



(19) 


3 Calculations 

3.1 Gaussian Charge Distribution 

One could assume, as a reasonable approximation to take both p^ and p-,, as Gaus¬ 
sian. That is 


PAr(u, CTat) — 




exp(- 


u 




P7r(v,a^) = 




<ixp(-^) (20) 


where \/<r‘^ >n = and \/< = -s/Sal- However, as a hrst approxi¬ 

mation we take the simpler case for which the meson is point-like, i.e., a.,r 0. Then 
P7r(v,cr7r) = 5 (v). The integral in Eq. ITSl reduces to 

bjv(u) 
r -|- u| 

and for p^ in Eq. [20] we have a rather simple result 



lh^ = lerf [ ^ 

+ u| r L\/2(TAr. 

Writing the Coulomb wave function appearing in Eq. [TK]as (pKci^) 
the equation for Rni{r) following from Eq. [THlis 



( 21 ) 

Rni{r)YiU0,<P) 


1 d ^^ dRnijr) 
dr dr 


—Rni{r) + 2p + «^erf 


_\/ 2 aN. 


Rni{r) = 0 


( 22 ) 
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For this approximation in which we take the meson radius to be zero and the 
nucleon rms to be \/3a^ we need to solve Eq. [22] and use the that we obtain 

to evaluate the contribution of the terms Eq. [11] where we can take v = 0 since 
p^(v) = (5(v). For this case, AE of Eq [19] becomes 




+ 


an 


1 


1 


-Q ■^°'N 


a 


N 


/ r/ \ 2 1 

2mambr ( V^uat ^ 


^ -A - 

■o-a ■ r X p)) 


r ■ (r ■ pjp - — 


-erf( 


) 


— (TAre 

n 


-e 


T \/2cTAr \/2(JAr 


(23) 


where the brackets mean that the values are to be calculated from the wave 
function obtained by solving Eq [22] 


contributions 

kaonic hydrogen 

kaonic helium^ 

kaonic carbon^^ 


2 

8616.7 

44739.9 

45876 


hnite size correction 

-1.24 

-295.3 

-61220 


total (Za)^ contribution 

8615.5 

44444.6 

393656 


Energy Contributions of order (Za)^ 

contributions 

kaonic hydrogen 

kaonic helium^ 

kaonic carbon^^ 

<p^> 

8 m 1 

-0.0302 

-6.0069 

-403.46 

<p*> 

8m7. 

-0.0044 

0.0326 

-0.03 


-0.0531 

-0.1699 

-0.50 


14.4623 

79.7487 

661.73 

■(■■■?)?> 

0.0011 

0.1326 

10.84 

< 2 mfm.r- ( ^ )erf( Jr ■ (r ■ p)p > 

0.0983 

0.7602 

0.28 

total (Za)^ contribution 

14.474 

74.4321 

268.31 


Energy Corrections order {Zay 


^The nuclear RMS used for the following tables was taken to be 1.2Y^3/5d^/^fm where A is the 
atomic number 
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4 Conclusion 


We have used the Bethe-Salpeter equation to calculate to (Za)^ electromagnetic 
contribution to the energy levels of kaonic hydrogen, kaonic helium (K —^ He) and 
kaonic carbon {K We have included hnite size effects by assuming a gaussian 

charge distribution in the heavier particle. To simplify the calculation we neglect the 
hnite size of kaon and calculate the much larger effect of the nuclear charge distribu¬ 
tion. As can be seen in the second table, these contributions are not negligible. All 
recoil contributions are obtained exactly by the use of the Bethe-Salpeter equation. 
The knowledge of the electromagnetic contribution to these energy levels is essen¬ 
tial to untangle the strong interaction from that of the electromagnetic contribution 

[1]. Furthermore the investigation of bound state systems such as these, because 
the interaction time is much longer than that in scattering experiments, allows for 
additional, precision tests of QCD. 
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